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1. Introduction 

The study of B physics is essential to determine the flavor structure of the Standard 
Model, through knowledge of the Cabibbo-Kobayashi-Maskawa (CKM) matrix de- 
scribing quark mixing and CP violation, which may be associated with the lack 
of symmetry between matter and anti-matter in the Universe. In fact, since the 
amount of baryons in the Universe predicted using the CKM mechanism is sev- 
eral orders of magnitude smaller than what is observed by astronomers, extensions 
of the Standard Model propose additional sources of CP violation, which must be 
tested against Standard-Model predictions. B mesons provide the ideal environ- 
ment for such tests.^ In particular, high-precision theoretical inputs are needed for 
hadronic matrix elements, which may be computed starting from the gauge theory 
itself using numerical simulations of lattice QCD. At present, however, it is not yet 
feasible to perform simulations on lattices that can simultaneously represent the 
two relevant scales of B physics: the low energy scale Aqcd, requiring large physical 
lattice size, and the high energy scale of the b-quark mass nib, requiring very small 
lattice spacing a. An approximate framework is therefore needed, but one should 
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strive to achieve sufficiently precise results, otherwise the task of overconstraining 
the parameters of the Standard Model is compromised. 

A promising such framework is to consider (lattice) heavy-quark effective theory 
(HQET), which allows for an elegant theoretical treatment, with the possibility 
of fully nonperturbative renormalization.'^'^ The approach is briefly described as 
follows. HQET provides a valid low-momentum description for systems with one 
heavy quark, with manifest heavy-quark symmetry in the limit mf, — !■ oo. The 
heavy-quark flavor and spin symmetries are broken at finite values of mi, respectively 
by kinetic and spin terms, with first-order corrections to the static Lagrangian 
parametrized by Wkin and Wgpin 



where 



£HQET ^ ^^{x)DoMx) - WkinOkin " Wspin Ospin , (1) 



Ckin = ^pf,{x)'D'^ Tph{x) , Ospin = ^|Jf,{x) a ■ B ijjhix) . (2) 



These 0(l/rn6) corrections are incorporated by an expansion of the statistical weight 
in l/jTib such that Okin, Cspin are treated as insertions into static correlation func- 
tions. This guarantees the existence of a continuum limit, with results that are 
independent of the regularization, provided that the renormalization be done non- 
perturbatively. 

As a consequence, expansions for masses and decay constants are given respec- 
tively by 



rUB = mbarc + E^'-' + C^kin^;^'" + ^.pin i^^^'" (3) 



and 



„ "mB ^HQET Stat ^1 , HQET SA , kin , spinN / a\ 

JB\I— = ^A P (1 + P + WkinP +UJspinP^ ), (4) 

where the parameters JTibaro and Z^'^^'^ are written as sums of a static and an 
0(l/mb) term (denoted respectively with the superscripts "stat" and "1/mf," be- 
low), and c^'^^^ is of order l/mt. Bare energies {E^^'^^, etc.) and matrix elements 
^pstat^ etc.) are computed in the numerical simulation. 

The divergences (with inverse powers of a) in the above parameters are can- 
celled through the nonperturbative renormalization, which is based on a matching 
of HQET parameters to QCD on lattices of small physical volume — where fine 
lattice spacings can be considered — and extrapolation to a large volume by the 
step-scaling method. Such an analysis has been recently completed for the quenched 
case.131 In particular, there are nonperturbative (quenched) determinations of the 
static coefficients Wbaro ^^'^ ^a'^^ fo'" HYPl and HYP2 static-quark actioni^l at 
the physical b-quark mass, and similarly for the 0{l/mb) parameters Wkin, Wspin, 

l/m, r.l/m, , HQET 
"^bare ' ^^^'^ ^A 

The newly determined HQET parameters are very precise (with errors of a 
couple of a percent in the static case) and show the expected behavior with a. 
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They are used in our calculations reported here, to perform the nonperturbative 
renormahzation of the (bare) observables computed in the simulation. Of course, 
in order to keep a high precision, also these bare quantities have to be accurately 
determined. This is accomplished by an efficient use of the generalized eigenvalue 
problem (GEVP) for extracting energy levels En and matrix elements, as described 
below. 

A significant source of systematic errors in the determination of energy levels in 
lattice simulations is the contamination from excited states in the time correlators 

oo 

C(t) = {0{t)Om = ^ |(n|d|0)pe-^"' (5) 

n=l 

of fields 0{t) with the quantum numbers of a given bound state. 

Instead of starting from simple local fields O and getting the (ground-state) en- 
ergy from an effective-mass plateau in C{t) as defined above, it is then advantageous 
to consider all-to-all propagators^ and to solve, instead, the GEVP 

Cit)v„it,to) = K{tM)C{h)vn{tM, (6) 
where t > to and C(t) is now a matrix of correlators, given by 

oo 

a,(i) = {O,{t)O,{0)) = *,i = i,---,^- (7) 

The chosen interpolators Oi are taken (hopefully) linearly independent, e.g. they 
may be built from the smeared quark fields using TV different smearing levels. The 
matrix elements are defined by 

EE (*„)i = {n\d,\0) , (m|n) = S^n ■ (8) 

One thus computes C'ij for the interpolator basis Oi from the numerical simu- 
lation, then gets effective energy levels and estimates for the matrix elements 
^ni from the solution A„(<,to) of the GEVP at large t. For the energies 

it is showrPthat E!^{t,tQ) converges exponentially as t — > oo (and fixed ip) to the 
true energy i?„. However, since the exponential falloff of higher contributions may 
be slow, it is also essential to study the convergence as a function of in order 
to achieve the required efficiency for the method. This has been recently done,^ by 
explicit application of (ordinary) perturbation theory to a hypothetical truncated 
problem where only N levels contribute. The solution in this case is exactly given by 
the true energies, and corrections due to the higher states are treated perturbatively. 
We get 



(^i^o) = En + en{t,to) 



(10) 
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for the energies and 

oo 

e-"'{Qf{t,to))^\0) = In) + 5] w(i,io)|n') (11) 

n' = l 

for the eigenstates of the Hamiltonian, which may be estimated through 

Qf{t,to)^Rnid,Vnit,to)), (12) 

In our analysis we see that, due to cancehations of i-independent terms in the 
effective energy, the first-order corrections in e„(t,io) are independent of to and 
very strongly suppressed at large t. We identify two regimes: 1) for to < V-^, the 
2nd-order corrections dominate and their exponential suppression is given by the 
smallest energy gap \Ern — En\ = AEm,n between level n and its neighboring levels 
m; and 2) for to > t/2, the Ist-order corrections dominate and the suppression 
is given by the large gap AE^+i^n- Amplitudes 7r„„'(t,to) get main contributions 
from the first-order corrections. For fixed t — to these are also suppressed with 
A£'7v+i,rf Clearly, the appearance of large energy gaps in the second regime im- 
proves convergence significantly. We therefore work with t, to combinations in this 
regime. 

A very important step of our approach is to realize that the same perturba- 
tive analysis may be applied to get the l/mt corrections in the HQET correlation 
functions mentioned previously 

a,it) = Qf*(t) + c.C^™^(t) + 0(0.2), (14) 

where the combined 0(l/mt,) corrections are symbolized by the expansion param- 
eter uj. Following the same procedure as above, we get similar exponential suppres- 
sions (with the static energy gaps) for static and 0(l/m6) terms in the effective 
theory. We arrive at 

(t,to) = Ef^^'-\t,to)+ujEf-'/"'^{t,to) + 0{uJ^) (15) 

with 

^eff,stat(^^ to) = + 13^1;^^ e-^-^"¥i." *-(-..., (16) 

^off,i/m,(^^^^) ^ ^1/™. + [^1/™. _ /3,f-ttA£;j^/;;''„]e-^^«" ."* + ... . (17) 

and similarly for matrix elements. Preliminary results of our application of the 
methods described in this section were presented recentlj^ and are summarized in 
the next section. A more detailed version of this study will be presented elsewhere. 1^ 



An (to + a, to) 
A„(to -I- 2a, to) 



(13) 
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2. Results 

We carried out a study of static-light Bs-mesons in quenched HQET with the non- 
perturbative parameters described in the previous section, employing the HYPl 
and HYP2 lattice actions for the static quark and an 0(a)-improved Wilson action 
for the strange quark in the simulations. The lattices considered were of the form 

X 2L with periodic boundary conditions. We took L « 1.5 fm and lattice spacings 
0.1 fm, 0.07 fm and 0.05 fm, corresponding respectively to (3 — 6.0219, 6.2885 and 
6.4956. We used all-to-all strange-quark propagators constructed from approximate 
low modes, with 100 configurations. Gauge links in interpolating fields were smeared 
with 3 iterations of (spatial) APE smearing, whereas Gaussian smearing (8 levels) 
was used for the strange-quark field. A simple 7075 structure in Dirac space was 
taken for all 8 interpolating fields. Also, the local field (no smearing) was included 
in order to compute the decay constant. 

The resulting (8 x 8) correlation matrix may be conveniently truncated to an N x 
N one and the GEVP solved for each N, so that results can be studied as a function 
of N. We then pick a basis from unprojected interpolators, by sampling the different 
smearing levels (from 1 to 7) as {1, 7},{1, 4, 7}, etc. We perform fits of the various 
energy levels and values of N to the behavior in Eq. and extract our results 
from the predicted plateaus. Next, we take the continuum limit, extrapolating our 
results to a — 0. We see that the correction to the ground-state energy due to terms 
of order l/m;,, which is positive for finite a, is quite small (consistent with zero) 
in the continuum limit. Our results for the pseudoscalar meson decay constant, 
both in the static limit and including 0(l/mb) corrections, are shown in terms 
of the combination (pf^Q^'^ ^ i^ps ^/m-Ps/Cps, where Cps(M/Aqcd) is a known 
matching function and ^^'-'^ de note s the renormalization-group-invariant matrix 
element of the static axial current .^Hl These two continuum extrapolations are shown 
in comparison with fully relativistic heavy-light (around charm-strange) datePl in 
Fig. [T] below. Note that, up to perturbative corrections of order in Cps, HQET 
predicts a behavior const. + O(l/rompg) in this graph. Surprisingly no l/(roTOps)^ 
terms are visible, even with our rather small errors. 



3. Conclusions 

The combined use of nonperturbatively determined HQET parameters (in action 
and currents) and efficient GEVP allows us to reach precisions of a few percent 
in matrix elements and of a few MeV in energy levels, even with only a moder- 
ate number of configurations. The method is robust with respect to the choice of 
interpolator basis. All parameters have been determined nonperturbatively and in 
particular power divergences are completely subtracted. We see that HQET plus 
0(l/mt,) corrections at the b-quark mass agrees well with an interpolation between 
the static point and the charm region, indicating that linearity in 1 /m extends even 
to the charm point. A corresponding study for Nf = 2 is in progress. 
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Fig. 1. Comparison of tlic continuum values for the pseudoscalar meson decay constant from Fig. 
4 to fully relativistic data in the charm region. The solid line is a linear interpolation between the 
static limit and the points around the charm-quark mass, which corresponds to l/rg mpg ~ 0.2. 
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